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The Brownian motion of a particle in a one-dimensional periodic potential subjected to 
a uniform external force F is studied. Using the formula for the diffusion coefficient D 
obtained by other authors and an alternative one derived from the Fokker-Planck equation 
in the present work, D is compared with the differential mobility fj, — dv/dF where v is the 
average velocity of the particle. Analytical and numerical calculations indicate that inequality 
D > fiksT, with /cb the Boltzmann constant and T the temperature, holds if the periodic 
potential is symmetric, while it is violated for asymmetric potentials when F is small but 
nonzero. 
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1. Introduction 

The response of a system in thermal equilibrium to an external disturbance has close 
relation to fluctuations produced spontaneously in the system in the absence the disturbance. 
This relation can be formulated as the fluctuation-dissipation theorem. 1,2 The Einstein rela- 
tion D = [loksT is a famous example, where D is the diffusion coefficient and hq is the mobility 
of a Brownian particle, k-Q is the Boltzmann constant, and T is the temperature. In this ex- 
ample, D measures the fluctuation of the particle position or velocity v and /j,q = limp^ov/F 
represents the response of the particle velocity to a small external force F. 

For systems far from thermal equilibrium, any particular relation between D and /xo 
is expected, because we do not know general laws like the fluctuation-dissipation theorem 
for such systems. However, recent investigations 3-5 into certain one-dimensional systems in 
nonequilibrium steady states suggest that inequality D > ^ik-^T with \i = dv/dF being the 
differential mobility may hold in these systems: numerical data show that D is greater than 
fJ-k^T for a Brownian particle moving in sinusoidal potentials, 3 for flushing ratchets 4 and for 
rocking ratchets. 5 Is there any rule that tells under what conditions inequality D > fik^T 
holds? Finding such a rule, if exists, would provide an important insight into understanding 
the behavior of nonequilibrium systems. 

The purpose of the present paper is to figure out whether inequality D > [ik^T holds 
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generally in the system of a Brownian particle moving in a one-dimensional periodic potential 
subjected to a uniform external force. This system is one of the simplest systems that exhibit 
nonequilibrium steady states, and convenient formulas for calculating D and (i are known. 3 ' 6 
From analytical and numerical investigations based on these formulas and the one we derive 
from the steady-state solution to the Fokker-Planck equation in the present work, we find 
that this inequality is likely to be valid for any symmetric potentials whereas it is violated for 
small external forces if the potential is asymmetric. 

2. Formulas 

We shall investigate the overdamped motion of a Brownian particle moving along the x 
axis under the influence of a periodic potential V(x) of period I and a uniform external force 
F. The total potential U(x) for the particle is given by 

U{x) = V(x) - Fx, V(x + l) = V{x). (1) 

In what follows periodic functions I±(x) = I±(x + /) defined by 

I±{x) = \ f e ±f>v{xWU(**v) dy (2) 
' Jo 

play important roles, where f3 = 1/ksT. The average of a periodic function f(x) of period / 
over the period will be denoted by (/): 

(f) = \ j\x)dx, f(x + l) = f(x). (3) 
The "normalized" functions 

J±(x) = I±(x)/(I±), (4) 

which satisfy (J±) = 1, are also of use. 

It was shown by Stratonovich 7 that the average velocity v of the particle can be calculated 
by the formula 3 ' 8 

v = D {l-e-P Fl )/l(I ± ), (5) 

where Do is the diffusion coefficient of a freely moving Brownian particle (V = F = 0) and it 
is related with the frictional coefficient £ of the particle through Do = ksT/^. It is noted that 
(7_l_) = The differential mobility fi = dv/dF can be calculated by differentiating eq. (5) 
with respect to F. The result can be expressed in a succinct form: 6 

fi = D (J + J_)/k B T. (6) 

The formula for D in the presence of both V(x) and F was derived recently by Reimann et 
al.: 3 

D = Do(J ± J+J-}. (7) 
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Note that (J+J_) is equal to (J+J 2 ). Reimann et al. 3 derived this formula by considering 
the moments of first passage time. Later, Hayashi and Sasa 6 obtained the same result by 
considering the system with an additional potential that varies much slowly than the original 
periodic potential V(x). 

If the periodic potential V(x) and the external force F are given, the diffusion coefficient 
D and the differential mobility /j, can be figured out by carrying out the two-dimensional 
integrals involved in eqs. (7) and (6); from the results we find whether or not D is larger than 
/ifceT. Nevertheless, an alternative formula may be useful in studying the sign of D — fik^T. 
From the steady-state solution to the Fokker-Planck equation we can derive (see the appendix) 
the formula 

D- nk B T = vl((J+-l)K-), (8) 

where periodic functions K±(x) of period I are defined by 

K±(x) = j £[J±(y)-l]dy. (9) 

Because the sign of v is the same as that of F as evident from eq. (5), formula (8) indicates 
that D > jikB,T if the sign of 

s = <(J+ - l)K.) (10) 

is the same as that of F. In analytic investigations, evaluation of eq. (10) is usually much 
easier than calculating eqs. (6) and (7) and then subtracting one from the other. By contrast, 
it is better to use eqs. (6) and (7) in numerical calculations, because the evaluation of the 
three-dimensional integral involved in eq. (10) is time consuming. 

3. Example 

In this section we present the numerical results for the diffusion coefficient D and the 
differential mobility fj, obtained from formulas (7) and (6), respectively, with a particular 
choice of potential: 

V{x) = A[sm(27rx/l) - Asin(47nz/Z)], (11) 

where A > and A are parameters. This potential is symmetric if A = or A = oo, and asym- 
metric otherwise. The potential height W, defined as the difference between the maximum 
and minimum values of V, is given by 

W = 2(1 - 2Ac)v / l - c 2 A, (12) 

where c is defined by 

c= (1- Vl + 32A 2 )/8A. (13) 

Note that this potential has a single minimum and a single maximum in each period if < 
|A| < 1/2, while it has an extra pair of local minimum and maximum ifl/2<|A|<oo. 
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Fig. 1. The difference D — [ik^T in units of D as a function of the external force F in the dimen- 
sionless form ((3FI) for the system with the potential given by eq. (11). The potential height is 
chosen as W — 5.0 k^T, and the results for different values of parameter A are shown. The inset 
presents the dependence of D/D and fikuT/Do on j3Fl in the case that (iW — 5.0 and A = 0.25. 



The inset of Fig. 1 shows the dependence of D and /u/cbT on the external field F in the 
case that (3W = 5.0 and A = 0.25. It appears that D is always larger than fiksT. However, 
closer inspection reveals that D is smaller than (ik^T in a certain range of F near F = 0: See 
Fig. 1, where the difference D — ^k^T is plotted against F in expanded scales for several values 
of A with (3W = 5.0; the results for negative values of A is obtained from the corresponding 
results for —A (which is now positive) by changing the sign of F, as the symmetry property 
indicates. In the case of symmetric potential (A = 0) we observe that D > fik^T (the equality 
holds when F = 0). 

For potentials with positive (negative) A, we find that D < fik^T in a range < F < Fq 
(Fq < F < 0) of F where the upper (lower) bound Fq depends on A, W and (3. Figure 2 
shows the dependence of (3FqI on f3W for several values of positive A. One sees that (3FqI is a 
monotonically increasing function of (3W. If A < 1/2 (the solid lines in Fig. 2), the value of 
(3FqI for a fixed (3W decreases with decreasing A and becomes zero as A = is approached. 
By contrast, (5FqI decreases with increasing A when A > 1 (the dashed lines in Fig. 2). These 
behaviors may be summarized that as the potential becomes symmetric, the range of F in 
which inequality D < [xk^T holds shrinks to zero. 

4. Conjectures 

We have calculated D and \x numerically for various periodic potentials V(x) in addition to 
the one described in the preceding section; some of the results will be presented in the following 
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Fig. 2. The upper bound Fq of the interval for external force in which inequality D < jik^T holds is 
plotted as a function of the potential height W for the system with the potential given by eq. (11). 
The results for different choices of parameter A are shown. 

section. We have also carried out analytic study on the sign of D — nk B T in several limiting 
cases, which will be discussed in the next section. From the results of these investigations, we 
have been lead to postulate the following conjectures. 

(i) Inequality D > [ik^T holds for arbitrary symmetric potentials V. 

(ii) Suppose that the potential is asymmetric and has a single minimum and a single maximum 
in each period. Let a be the distance from a minimum to the adjacent maximum on the 
right (see Fig. 3). Then, we have D < ^k B T for < F < F (F < F < 0) and D > ^k B T 
outside this interval of F if a > 1/2 (a < 1/2), where Fq is a positive (negative) constant 
that depends on potential V(x) and temperature T. 

Note that in the example considered in the preceding section distance a is given by 



where c is defined by eq. (13), and hence condition a > 1/2 corresponds to A > 0. Therefore, 
the results shown in Fig. 1 are consistent with these conjectures. 

5. Evidence 

The conjectures stated in the preceding section are based on the analyses presented in 
this section. We first describe the analytical investigations, in several limiting cases, into the 
sing of D — fiksT using formula (8). Then, considering the results of these investigations and 
supplementary numerical calculations, we will argue for the validity of the conjectures. 



a = (I /it) arccos c, 



(14) 
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Fig. 3. Two examples of periodic potential V(x) of period I that has a single minimum and a single 
maximum in a period are schematically shown. The one represented by the solid line has rounded 
peaks at its maxima and rounded valleys at its minima, while the one represented by the dashed 
line has cusps at its maxima. The location of a minimum is xq, and the the distance from this 
minimum to the adjacent maximum on the right is a. The potential height W is defined as the 
difference between the maximum and minimum values (V m&K and Vmi n ) of V. 



5. 1 Small external force 

The first limiting case we study is the case of small external force represented by condition 
(3\F\l <C 1. In this case the factor s defined by eq. (10) may be expanded in powers of (3FI as 

s = so + SiPFl + s 2 {l3Flf + . . . . (15) 

In order to express the expansion coefficients sq, s\ and so on concisely, we introduce periodic 
functions (of period I) ip±(x) and X±( x ) by 

^ ± ( x ) = e ±/3V »/(e ±/3y > (16) 

and 

X±(x) = ]£[^±(y)-l]dy. (17) 
It is not difficult to see that periodic functions I±(x) defined by eq. (2) can be expressed as 
I±(x) = (e^V^jl - iiFl [1/2 T Xt (x) ± <X T >] 

+ 0[{f3Fi) 2 ]}. (18) 
Therefore, the normalized functions J±{x) defined by eq. (4) are given by 

J±(x) = ip±(x){l ± m lx T (x) - (V±Xt)] 

+ 0[(f3Fl) 2 ]}, (19) 
from which the following expression for K±{x) defined by eq. (9) is obtained: 

K±(x) = X ±(x) ± pF ! ^±{y)[xT(y) ~ <V>±X T >] dy 
J o 

+ 0[{I3FI) 2 }. (20) 
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Substituting eqs. (19) and (20) into eq. (10), one finds 



«o = ((V>+ - 1)X-) 



(21) 



and 



si = (tp+x 2 -) - (4> + x-) 2 + W-xl) - <V>-x+> 2 - 



(22) 



It is worth noting that coefficient si cannot be negative: 



si > 



(23) 



where the equality holds only in the trivial case of a constant potential V. This property 
comes from the Schwarz inequality 



and identity (tp±) = 1 resulting from the definition (16) of ip±; the equality in eq. (24) holds 
if and only if x± is a constant (i.e., V is a constant). 

By contrast, the leading term so in expansion (15) can be positive or negative. However, in 
the case of symmetric potential, i.e., if there exists a constant a such that V{x) = V(2a — x) 
holds for any x, we have sq = 0. The reason is the following: for such a symmetric potential, 
4>±(%) is symmetric and Xt( x ) i s antisymmetric about x = a, hence we obtain (iP±Xt) = and 
(Xqr) = 0, which imply so = 0. This fact and inequality (23) indicate that inequality D > fiksT 
holds for any symmetric potentials as long as F is small, which supports conjecture (i) stated 
in the preceding section. 

In the case of asymmetric potential, it is expected that sq ^ 0. Then, what property of V 
determines the sign of sq (i.e., the sign of D — fik^T for small F)? It seems difficult to answer 
this question for arbitrary potentials. However, if we restrict our attention to a certain class 
of potentials, we can, at least partly, answer the question. Let us consider a potential that has 
only one minimum and one maximum in a period, as shown in Fig. 3. Let xo be the location 
of a minimum, a be the distance from a minimum to the adjacent maximum on the right, W 
be the potential height defined as the difference between the maximum and minimum values 
of V. The potential V may have a rounded peak at its maximum and a rounded valley at its 
minimum as shown by the solid line in Fig. 3. It may have a cusp at its maximum (the dashed 
line in Fig. 3), or at its minimum, or at both. We shall analyze the sign of sq in the limiting 
cases of large potential height (f3W S> 1) and small potential height {[5W <C 1). 

Let us consider the case of large potential height, (3W S> 1. In order to make the analysis 
simple, we assume that the origin of the x axis is chosen in such a way that condition < xq < 
xo + a < I is satisfied. In evaluating sq given by eq. (21), it is noted that function tp+(x) has 
a sharp peak at x = xq + a and vanishes rapidly as one moves away from the peak. Therefore 
so can be approximated by 



(24) 



so - (ip+)X-{xo + a) - (x_), 



(25) 
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since x~( x ) does not vary rapidly in the vicinity of x = xq + a as we shall see in a moment. 
Function x~( x ) defined by eq. (17) is the sum of 

v{x) = \ J%-(y)dy (26) 

and —x/l. Since the integrand ~4>-(y) in eq. (26) is practically zero except a narrow region 
around the sharp peak at y = xo, function a(x) behaves like a step function: as x is increased 
from zero to I, a(x) increases rapidly from zero to unity around x = xq. Therefore, X-( x ) is 
well approximated by X-( x ) — 1 — x /l near x = xo + a and it does not change rapidly in the 
vicinity of x = xq + a. We also find that (%_) = 1/2 — xq/1 if the small correction of order 
1//3W is neglected. From these arguments and identity = 1 we obtain 

*4-f (27) 

This expression for sq reveals that the sign of so is determined by whether the location of 
the top of the potential hill between a pair of neighboring valleys is closer to the left valley 
(a /I < 1/2) or the right one (a /I > 1/2), which supports conjecture (ii) in the preceding 
section. 

Now we turn our attention to the case of small potential height, [3W <C 1. It will be 
assumed that (an arbitrary constant is added to V such that) the maximum and minimum 
values of V are of order W. Then condition (3W <C 1 implies \(3V\ <C 1. Expanding ip±(x) 
and x±( x ) defined by eqs. (16) and (17) in powers of (3V, and then substituting them into 
eq. (21), we obtain 

so = -P 3 ((V 2 - (V 2 ))V) + 0[(f3V)% (28) 

where periodic function V(x) of period I is defined by 

V(x) = jJ*[V(y)-{V)]dy. (29) 

Unlike the case of (3W » 1, we have not been able to relate the sign of so approximated 
by eq. (28) to that of 1/2 — a for general asymmetric potentials. Here we investigate the sign 
of so for three examples of potential V(x). The first example is the one considered in Sec. 3, 
see eq. (11). The second example is a piecewise-cubic function given by 

V{x) = A{{x/lf + X{x/l)[l-A{x/lf]} (|x|<//2), (30) 

where A and A are parameters; V(x) outside the range \x\ < 1/2 is defined such that it is a 
periodic function of period I. We shall assume that A > and |A| < 1/2. Then V(x) has a 
cusp at its maximum, as the one represented by the dashed line in Fig. 3 does. The distance 
a from a minimum to the adjacent maximum on the right is given by 

a = l(- + A ^ , (31) 

V2 1 + V1 + 12AV V ^ 
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Table I. Approximate expression for s given as the leading term in cq. (28) obtained for three 
examples of V(x). 





Example 1 


Example 2 


Example 3 


V(x) 


eq. (11) 


eq. (30) 


eq. (33) 


so 


3A(/3A) 3 
16tt 


A(/3A) 3 (i l x 2 n 
1575 V12 11 > 


360 V 1 l ) 



and the potential height W by 



W = A 



2 ^ 1 + 12A 2 



_9 18(1 + V1 + 12A 2 ). 
The third example is a piecewise-linear (sawtooth) potential 



(32) 



)Wx/a < x < a 

(33) 
W{l-x)/(l-a) a<x<l, 

where W and a are positive parameters with restriction < a < I; again, V(x) outside the 
range < x < I is defined such that it is a periodic function of period I. Parameter W 
represents the potential height, and parameter a corresponds to the distance from a minimum 
of V(x) and the adjacent maximum on the right. 

For each example, the leading term of so given in eq. (28) has been calculated. The results 
are summarized in Table I. In all the three examples the sign of so is the same as that of 1/2 — a 
(remember that //2>aifA<0in the first two examples). This observation is consistent 
with conjecture (ii). 

It is interesting to note that sq is of higher order in (3W than 

S1 = 2(3 2 {(V 2 }-(V} 2 )+01(PV) 3 } (34) 

in the case of small potential height. This fact implies that s change its sign at small (3FI 
when the latter is varied. Let Fq be the value of F at which s changes its sign, then one finds 
from eq. (15) that 

mi ~- s ^~ m " ~ {y2))v) (35) 

W - ai - 2«V 2 > - (V) 2 ) ' (35) 
which is of order (3W '. For the first example considered above, we obtain 

3-7rA 

m * wrm/ A - (36> 

This relation and eq. (12) explain the behavior of the graphs in Fig. 2 near the origin. 
5. 2 Large external force 

If the external force F is large enough (f3\F\l S> 1), the dominant contribution to the 
integral in eq. (2) defining I±(x) comes from the narrow region near y = (if F > 0) or y = I 
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(if F < 0). Therefore, I±(x) may be expanded as 

/ ± ( x )~£__y [/ lT (0)+^(0) + ...]e-^dy (37) 

if F > 0, where hzp(y) = e^^^^^, and h'^y) is the derivative of hzp(y). A similar expansion 
in the case of .F < can be made. Using these expansions, periodic functions I±(x), J±(x), 
and K±(x) are expressed as the power series in 1/F. Substitution of J+(x) and K-(x) thus 
obtained into eq. (10) yields 



1 

s = 



f3Fl 



2{(V'Y) , 5((n j ) , ( 1 



(38) 



^2 1 j?3 I J?^ 

where V'(x) is the derivative of potential V(x). This expression is valid both for F > and 
for F < 0. The leading term of s given by eq. (38) has the same sign as that of F and hence 
inequality D > fik^T holds if \F\ is large enough. 
5.3 Small potential height 

The last limiting case we study is the limit of small potential height; the strength of the 
external force F can be arbitrary. In this case we find it convenient to express the potential 
V{x) in the Fourier series as 

°° 2tt 
V(x)= Y, Vne iknX , k n = —n. (39) 

n=— oo 

In the integrand of eq. (2), factor e^@ vlyX ^ y ^ is expanded in powers of (3V(x =F y) an d then 
eq. (39) is substituted to carry out the integral. Once I±(x) are obtained in this way, it is 
straightforward to calculate J±(x) and K±{x). Substituting the resulting expressions for J+(x) 
and K-(x) into eq. (10), we have 

The sign of the leading term in this expression for s is the same as that of F, and hence 
inequality D > /ifesT holds if (3W is small enough. 

It is noted that in the limit of small (3FI the leading term in eq. (40) approaches to 

00 \RV I 2 

^ Fl E TfT^ = 2 ^(( v2 ) " (V) 2 )/5^- (41) 

n=l { n) 

This expression agrees with the leading term of eq. (34) multiplied by (3FI. This is expected 
from the consistency of the analysis. Similarly, the term of order ((3V) 3 in eq. (40) should 
converge to the first term of sq given in eq. (27) in the limit F — > 0, which we have not 
checked. In the opposite limit, (3\F\l » 1, the leading term in eq. (40) converges to the 
leading term in eq. (38), because J2n^n\^n\ 2 = ((U') 2 )- 
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5.4 Symmetric potentials 

Here, we consider the case of symmetric potential and argue for the validity of conjecture 
(i). In this case, s is an odd function of F (D — /x/cbT is an even function of F), and therefore 
we need to examine the sign of s only for F > 0. Remember that s > is equivalent to 
D > fikftT when F > 0. It has been shown that 

s ~ Sl f3Fl (42) 

with si > for small 0FI (§5.1) and s ~ 2((V') 2 ) / pF 3 l for large /3FZ (§5.2). Hence, it is 
concluded that inequality D > [ik^T holds (the equality holds when F = 0) in these two 
extremes. Furthermore, this inequality has been found to be valid in the entire range of F if 
the potential height is small compared to the temperature (§5.3). 

In order to assert the validity of conjecture (i), we have to demonstrate that s > for 
intermediate values of (3FI when f3W is not small. For this purpose, numerical calculations of 
s = (D — fik^,T)/vl are carried out using formula 

s = <4M -(/+)(/+/-) (43) 

<J+) 2 (l-e-^) 1 j 

obtained from (5), (6) and (7); as remarked earlier, this method of evaluating s is more 
convenient for numerical calculations than using formula (10). Symmetric potentials V(x) of 
the following type are examined: 

N 

V{x) = A c n cos(2n7rx/Z), (44) 

n=l 

where is a positive integer, c n are arbitrary coefficients, and the overall factor A is deter- 
mined such that the potential height is W for given values of W and c\, C2, . . . , cat. 

Figure 4 shows the numerical results for a potential with an arbitrarily chosen set of co- 
efficients {c n } in the case of N = 3. Here, s is plotted as a function of (3FI for several values 
of (3W . As F is increased from zero, s starts to increase linearly in F as eq. (42) predicts and 
it continues to increase until it reaches a maximum value, and then decreases monotonically. 
Qualitatively the same behavior of s are observed for other potentials corresponding to differ- 
ent sets of {c n } with = 3 or = 5 (data not shown), which strongly suggests the validity 
of conjecture (i). 

In Fig. 4, the analytic results, the leading terms in eqs. (40) and (38), are also plotted. It 
is remarkable that the approximate expression (40), which is valid for small (3W, agrees quite 
well with the numerical results for [3W as large as f3W ~ 1. For [3W larger than about unity, 
the dependence of s on F is well approximated by the leading term of eq. (38) if F is larger 
than a few to several times the maximum slope V^ ax = m&x. x {V'(x)} of potential V(x); in 
the example shown in Fig. 4, V^ax ^ 8AW/1. 

In addition to the numerical analysis concerning the dependence of s on F, shown in 
Fig. 4, for more than ten different potentials, we have carried out more extensive search for 
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Fig. 4. Dependence of s on (3FI numerically obtained for a symmetric potential (44) with N = 3 
and ci = 0.265947, c 2 = 0.823433, c 3 = -0.522984; the inset depicts the potential function. 
The results for different choices of potential height W are shown. The dashed line represents the 
analytic expression (40) valid for small (3W, and the dash-dotted lines indicate the leading term 
in eq. (38) for large f3Fl. It is remarked that the graphs of 10 times s instead of s are plotted for 
f3W = 1.0. 



possibility of negative s. Symmetric potentials expressed by eq. (44) with N = 3 and those 
with N = 5 are studied. For a given N, every coefficient c n (n = 1, 2, . . . , N) is chosen from a 
random number uniformly distributed in interval (—1, 1). For each set {c n } of the coefficients, 
the potential height W is chosen from a uniform random number in interval (0, W max ) where 
Wmax is set to be 20k^T; and for a given W, the external force F is chosen from a uniform 
random number in (0, F max (W)) where -F max is set to be 2V^ ax . We have examined 1000 sets 
of {c n } and 300 sets of {W, F} for each set of {c n } in the case of N = 3, and 2000 sets of {c n } 
and 100 sets of {W, F} in the case of N = 5. In the data of these 5 x 10 5 samples we have not 
detected any instance in which s < 0. 

All these analytical and numerical investigations firmly indicate that the statement of 
conjecture (i) should be true. 
5.5 Asymmetric potentials 

Now we discuss conjecture (ii) associated with asymmetric potentials. If the potential 
height is small {(3W <C 1), the analyses of §5.1 and §5.3 show that inequality D > \xk^T holds 
for almost entire range of F except a small interval of order W/l. This interval is given by 
< F < Fq or Fq < F < depending on the sign of Fq given by eq. (36). 

If the potential height is not small, we do not have enough evidence to support conjecture 
(ii). It is true that D — fiksT change its sign at F = when F is varied (§5.1) and that it is 
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positive for large enough \F\ (§5.2). Furthermore, it is shown (§5.1) that in the case of large 
potential height (f3W > 1) we have D < fik B T for F > (F < 0) if a > 1/2 (a < 1/2) and 
\F\ is small. These results are consistent with conjecture (ii), but we are not certain, from 
the analytical study given above, whether there is only one interval on the F axis (as the 
conjecture states) where inequality D > nk B T is not satisfied. The numerical investigation 
presented in §3 for potential given by eq. (11) and a similar one (data not shown) for the 
piecewise-linear potential (33) support the validity of conjecture (ii). 

6. Conclusion 

We have postulated two conjectures (§4) concerning the diffusion coefficient and the dif- 
ferential mobility of a Brownian particle moving in a one-dimensional periodic potential under 
the influence of a uniform external force. We are quite certain about the validity of conjecture 
(i) associated with symmetric potentials (§5.4). It should be possible to prove it mathemati- 
cally, although we have not yet succeeded. Conjecture (ii) related with asymmetric potentials 
is partly speculative (§5.5). 

The ratio 6 = D/jik^ may interpreted as an effective temperature 5 ' 6 of the system in 
nonequilibrium steady state. Then, our conjectures imply that the effective temperature is 
higher than the temperature of the heat bath if the potential is symmetric or if the external 
force is not too small in the case of asymmetric potential. 

Very recently Hayashi and Sasa 10 have reported an alternative inequality associated with 
the diffusion coefficient and the differential mobility. They have proved that inequality D/Dq> 
{lik^T / Dq) 2 holds in general for the system considered in the present work. 
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Appendix: Derivation of eq. (8) 

Our derivation of formula (8) is based on a prescription to calculate the diffusion coefficient 
from the solution to the Fokker-Planck equation. 4 ' 8 ' 9 Let P(x) be the probability distribution 
function of the particle in the steady state. It satisfies the Fokker-Planck equation 



for the steady state. We assume that P{x) is periodic [P(x + /) = P (x)] and normalized such 
that Jq P(x) dx = 1. Such a solution is found to be given by 




(A-l) 



P{x) = J-{x)/l. 



(A-2) 



The average velocity v can be calculated from P(x) as 




(A-3) 
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and this leads to formula (5). Note that the right-hand side in eq. (A-3) is independent of x 
due to the Fokker-Planck equation (A-l). In order to calculate the diffusion coefficient D, we 
need to solve the differential equation 

^(e + ^)«« = (' + *s)™-t <a ' 4) 

for Q(x), where P(x) is the probability distribution function given by eq. (A-2). The diffusion 
coefficient is calculated from a periodic solution Q(x) = Q(x + I) to eq. (A-4) as 

D = D - (PD ^- + v^j Q(x) dx. (A-5) 

Any periodic solution Q yields the same result for D. Festa and d'Agliano 9 solved eq. (A-4) 
in the case of no external force (F = 0), and obtained a formula for D, which is similar to 
eq. (7) but much simpler. Here, we solve eq. (A-4) in the case of nonzero external force, and 
derive eq. (8). 

Integrating eq. (A-4) once, we have 

jL + ^)g (l) = ,(,), (A.6) 

where q(x) is given by 

q(x) = J-(x)/l + vK_(x)/D Q . (A-7) 

Here, K-(x) is defined in eq. (9). We have chosen the integration constant arbitrarily to get 
q{x) in eq. (A-6), since any periodic solution Q(x) is acceptable as remarked above. Integrating 
eq. (A-6), we arrive at 

„-/3U(x) rl 

= ' i-e-m J ^ U{x+y) q(x + y) dy (A-8) 

after some manipulations. This time, the integration constant has been determined such that 
Q(x) is periodic. 

Now we substitute eq. (A-8) into eq. (A-5) to study the diffusion coefficient. Making use 
of eq. (A-6) and the periodicity of Q(x), we rewrite eq. (A-5) as 

D = D — D [ q(x) dx-v [ Q(x) dx. (A-9) 
Jo Jo 

The second term, without the minus sign, on the right-hand side in this equation reads 

D / q(x)dx = D + vl{K-), (A- 10) 

Jo 

according to the definitions of q(x) and J_(x). Insertion of eq. (A-8) into the third term in 
eq.(A-9) yields the integral 

fdx e -P u W tdy e^ U( - x+y) q(x + y) = I f I + {x)q{x) dx, (A-ll) 
Jo Jo Jo 
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where the right-hand side is obtained by interchanging the order of integral and by using 
the fact that U(x + y) — U(x) and q(x) are periodic functions of x. From this identity and 
eqs. (A-8) and (5) we have 



' / Q(x)dx = — Do / J + (x)q(x) dx 
Jo Jo 



= -nk B T - vl{J+K-), (A- 12) 

where the second equality is due to eqs. (A-7) and (6). Substitution of eqs. (A-10) and (A-12) 
into eq. (A-9) gives eq. (8). 

The equivalence between formula (8) for the diffusion coefficient and the one, eq. (7), 
obtained by other authors can been shown as follows. Since it can be seen by integration by 
parts that ((J + — l)K + ) = 0, eq. (8) may be written as 

D = fiksT + vl((J + - l)(K- - K + )). (A-13) 

Now, it is not difficult to see from the definitions of J±(x) that 

l WM= "-^-W . (A,4) 
This relation and the definition (9) of K±(x) lead to 

icy - k +(x ) = tMfeMzWTO. (A . 15) 

Substituting this equation into eq. (8) and using formula (5) for v, we find 

D = f ik B T + Do(JlJ-)-D (J + J-). (A-16) 

Here, the first and the third terms on the right-hand side cancel out due to eq. (6). Therefore 
eq. (A-16) is identical to formula (7), and the equivalence between eqs. (8) and (7) has been 
proved. 
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